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Abstract. On January 6 th 2004, the IBAS burst alert system triggered the 8 th gamma-ray burst (GRB) to be located by the 
INTEGRAL satellite. The position was determined and publicly distributed within 12 s, prompting ESA's XMM-Newton to ex- 
ecute a ToO observation just 5 hours later, during which an X-ray afterglow was detected. The GRB had a duration ~ 52 s with 
two distinct pulses separated by ~ 42 s. Here we present the results of imaging and spectral analyses of the prompt emission from 
INTEGRAL data and the X-ray afterglow from XMM-Newton data. The y-ray spectrum is consistent with a single power-law 
of photon index -1.72 ± 0.15. The fluence (20 - 200 keV) was 8.2 x 10~ 7 erg cm -2 . The X-ray afterglow (F,,(r) oc v^T 6 ) was ex- 
tremely hard with fix = 0.47±0.01 and 5 = 1.46 ±0.04. The 2- lOkeVflux 1 1 hours after the burst was 1.1 x 10~ 12 erg cm" 2 s -1 . 
The time profile of the GRB is consistent with the observed trends from previous analysis of BATSE GRBs. We find that the 
X-ray data are not well-fit by either a simple spherical fireball or by a speading jet, expanding into a homogeneous medium or 
a wind environment. Based on previously determined correlations between GRB spectra and redshift, we estimate a redshift of 
~ 0.9^Q 4 (1 o") and a lower limit on the isotropic radiated energy of ~ 5 x 10 51 erg in this burst. 

Key words, gamma-rays: bursts - gamma-rays: observations - X-rays: observations 



03 1. Introduction l2003h . coupled with two monitors, one in the X-ray band and 

>the other working at optical wavelengths. Although not built as 
Gamma-ray bursts are an amazingly energetic phenomenon, ^r>o • * a ■ ■ nnr^Dif u u * i * 

•r-j , ■> si a GRB-onented mission, INTEGRAL has a burst alert system, 

^ . capable of a jetted output of order 10 51 erg in a few seconds. ffiAS (Mere2hetti et al 2003) ffiAS carries out rapid locali . 

H ; Although first detected in the late 1960s, significant progress sations for QRBs incident on the IBIS dete ctor with a preci- 

^ , has mostly been achieved in the last dozen years. That GRBs sion of a few (Mereg hetti et al. 2004b). The pub- 

are extragalactic in origin was suggested by the isotropic dis- ,• . .• f . u ,■ . U1 ,.. , 

& & && 1 F he distribution of these co-ordinates enables multi-wavelength 

tribution of GRBs observed by BATSE on board the Compton , e e . , . , ■ nrTO^oo a * 

j — ' — — it -| | — «- — | searches for afterglows at lower energies. INTEGRAL data 

Gamma-Ray Observatory ^ Meegan et alJ | l992t | Fishman et al. | Qn ^ emission ^ combination with the early multi . 

1994). The discovery by BeppoSAX of afterglows in the X- wavelength studies offer the best currently ava0able probe of 

ray dCosta et al. 1997) and subsequent discoveries at optical .. , 

j ' 1 1 II i M j ; 1 1 3 F these energetic phenomena. 

( van Paradiis et al. 1997) and radio (Frail et al. 1997) wave- 
lengths have led to redshift measurements jMetzger et"a"fl 

1 19971) for ~ 40 bursts ranging from 0. 105 to 4.5. 2. INTEGRAL Observations and Results 



ESAs International Gamma-Ray Astrophysics Laboratory GRB 040106 wafj detec ted by IBAS at 17:55: 1 1 UTC on 
INTEGRAL J Winkler et alj | 2003|) , launched in October 2002, January6 th 2004 , Mereg hetti et alJ EooS) with a signal to 
is composed ||0 f two main telescopes, an imager IBIS noifje ratio of ? ? jn ^ 15 tQ 200keV band for the inter . 

val from 17:55: 10UTC to 17:55:65 UTC. The burst was ob- 
,h 



JUbertini et alJl2003l) . and a spectrometer SPI ( IVedrenne et alJ 



* Based on observations with INTEGRAL, an ESA project with served at ™0 = H" 52™ 17. I s , 6 J2(m = -46° 47' 15", at an 

instruments and science data centre funded by ESA member states off-axis angle of 10.5°, lying in the partially coded FoV of 

(especially the PI countries: Denmark, France, Germany, Italy, IBIS and of SPI, but outside the FoV of the two monitoring 

Switzerland, Spain), Czech Republic and Poland, and with the par- instruments, JEM-X and the OMC. The IBAS alert was au- 

ticipation of Russia and the USA. tomatically distributed approximately 12 s after the burst start 
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time with a positional uncertainty of only 3'. Due to the weak- 
ness of the burst, data in the energy range 20-60keV from 
the two time intervals around the prominent peaks of emis- 
sion were combined to enable SPI to determine a position for 
the GRB. The position for GRB 040106 extracted from the 
SPI data is 0,2000 = ll*52 m 51', (5/ 2 ooo = -46° 47' 13"(S/N = 
6.8), which is 5.7' away from the IBIS location, consistent 
with the 10' location ac curacy of SPI for a source with S/N 
~ 10 JPubafh et alJ2004 . At 23:05 UT XMM-Newton began a 
45 ks exposure ToO observation. A bright source was visible 
at Qr/2000 = 1 \ h 52'" \2 A\ S.nnm = -46° 4 7' 15.9" in the 30 ks 
Quick-L ook- Analysis (Ehle et alJl2004T) . 0.9' from the IBAS 
position jTedds & Wa tson 2004). 

2.1 . Light Curves 

This GRB falls into the class of long bursts with a duration 
of ~ 52 s and two prominent pulses with a peak to peak sep- 
aration of 42 s and a long quiescent interval of ~24s be- 
tween the pulses. Light curves a re availa ble from SPI and 
IBIS. The ISGRI detector of IBIS jLebrun et alJl2003l) . an ar- 
ray of 128 x 128 CdTe crystals sensitive to lower energy y— 
rays, was used to produce the light curves of GRB 040106 
in the ranges 15 -40keV and 40-200 keV (Fig.Q}. The burst 
light curves were extracted using the pixels which had at least 
half of their surface illuminated by the GRB. The spectral lag 
could not be determined with any precision due to the weak- 
ness of the burst. The best constraint was obtained using a 
cross-correlation function on the background-subtracted IBIS 
data in the energy ranges 15-40keV and 40-200keV and a 
time-binning of 0.2 s from 10-65 s after 17:55:00 UT. The mea- 
sured lag was ~ + 1 s which is about what can be estimated by 
eye and does not allow us to constrain the GRB to a location 
on the lag-luminosity diagram, since lags are expressed on a 
logarithmic scale. 

The determination of SPI light curves with a short time 
binning is only possible by using the 1 s count rates of the 
19 germanium detectors, which are usually used for scientific 
house-keeping purposes. These values reflect the count rates 
of each detector in the broad SPI energy band from 20 keV - 
8MeV. The SPI light curve in Fig. was generated by sum- 
ming the background-subtracted count rates. The background 
(~ 50 counts s _1 detector -1 ) was determined from a 40 minute 
period before the burst trigger and was subtracted for each de- 
tector individually. 

A hardness ratio was derived from the ISGRI data by com- 
paring the count rate in the different energy bands such that 
HR = (H-S)/(H + S). The time bins were chosen to ensure at 
least 100 counts per bin in the sum of the two bands. The evolu- 
tion of the hardness ratio is shown in the bottom panel of Fig. ^ 
Although visual inspection indicates a slight hardening in the 
second pulse, this is not statistically significant (2 cr level). 

The I BIS light c urve w as denoised using a wavelet anal- 
ysis (Ouill iganetal.ll2002l) and the temporal properties of the 
two pulses (rise time, fall time, FWHM) were extracted and are 
shown in Table 
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Fig.l. Light curves of GRB 040106 beginning at 17:55 UTC 
extracted from (top to bottom) IBIS/ISGRI data in the range 
15-40keV, IBIS/ISGRI data in the range 40-200keV, SPI 
scientific housekeeping data in the range 20 ke V - 8 MeV and 
(bottom panel) the hardness ratio evolution derived from 
IBIS/ISGRI count rates. 



Table 1. Temporal properties of the pulses in GRB 040106. 



Parameters 


Rise Time, t r 


Fall Time, tf 


FWHM 




(s) 


(s) 


(s) 


1 st Pulse 


4 


10 


5 


2 nd Pulse 


10 


6 


5 



2.2. Spectral Analysis 

IBIS is a coded mask instrument and the photons of a sin- 
gle point source are spread over all the individual detectors. 
Spectral extraction is possible using specifically designed soft- 
ware which consists of modelling the illuminated mask by a 
point source of unitary flux placed at the sky coordinates of 
the GRB. The model is then fit to the detected shadowgram 
in each energy channel to obtain the rate and error for each 
channel. ISGRI single events were used to derive photon in- 
dices and fluxes for two time intervals, of 9 s and 8 s dura- 
tion, around the peaks of emission as shown in the column 
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Table 2. Spectral analysis results for GRB 040106 with SPI and 
IBIS/ISGRI for two intervals around the prominent peaks of 
emission. Errors quoted are for 1 parameter of interest at lcr 
confidence level. Fluxes are quoted for the energy range 20 - 
200 keV. Within the large errors on the spectral fits, results are 
consistent between instruments. 



Inten'al 


Parameter 


SPI 


IBIS / ISGRI 


1 st pulse 

(9 s) 


Photon Index 
Flux (erg cm -2 s~ 


-1 47 +0 - 52 
l ) 5.9 x 10- 8 


-1.70 ±0.12 
2.5 x 10~ 8 


2 nd pulse 
(8 s) 


Photon Index 
Flux (erg cm -2 s~ 


i Q9+O.31 
1-J -0.34 

l ) 6.1 x 10- 8 


-1.41 ±0.1 

3.2 x 10~ 8 



'IBIS/ISGRI' of Table El Spectral fitting was carried out using 
two methods. Firstly, a time averaged spectrum was derived 
for GRB 040106 by comparing its count rate in different en- 
ergy bins to the corresponding values obtained from the Crab 
Nebula observed at a similar position in the FoV. This method 
yields a best-fit single power-law model with photon index 
-1.32 + 0.28 and a fluence of 1.1 x 10~ 6 ergcirT 2 . Secondly, 
the effect of the IBIS mask support structure (the so-called 
nomex), which absorbs low-energy (<50keV) photons differ- 
ently at different off-axis angles has been taken into account us- 
ing Offline Software Analysis (OSA) version 4.0 off-axis cor- 
rection matrices and the extracted spectra have been modified 
accordingly. The improved low-energy response yields a best- 
fit single power-law model with photon index of -1.72 + 0.15 
and a fluence of 8.2 x 10~ 7 erg crrT 2 . The spectral fit results for 
the 2 pulses using the nomex correction are given in Table |3 
Due to the improved low energy response, these are the spec- 
tral indices used in subsequent discussion. 

The spectral evolution of GRB 040106 was also investi- 
gated with SPI. Spectra were extracted for two time inter- 
vals around the peaks of emission. The background was de- 
termined from a 40 minute period of SPI data before the burst 
trigger. The first interval is 7 s long starting at the very be- 
ginning of the burst at 17:55: 11UTC, and the second be- 
gins 34 s later and lasts for 12 s. Single events detected by 
SPI, corrected for intrinsic deadtimes and telemetry gaps, were 
binned into 5 equally spaced logarithmic energy bins in the 
20keV to 500 MeV range in each of the chosen intervals. 
Version 3 of the OSA fro m the INTEGRAL Science Data Centre 
ll Courvoisier et al.ll2003h a nd the software package SPIROS 6 
( Ski nner & Co nnell 2003) were used for SPI spectral extrac- 
tion, while XSPEC 1 1.2 was used for model fitting. The results 
are shown in the column labelled 'SPF in Table |2 The same 
analysis was conducted for multiple events incident on the SPI 
detectors, but yielded no significant improvement to the fit. 

3. XMM-Newton 

XMM-Newton observed the position of GRB 040106 for about 
45 ks, starting only ~ 20 ks after the burst. For the MOS1 cam- 
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Time since burst (s) 

Fig. 2. Background subtracted EPIC light curve of the X-ray 
afterglow of GRB 040106 in the 0.4 - lOkeV energy band. The 
dashed line shows the best-fit power-law decay 5=1 .46. 



era jTurner et alJ200ll) of the EPIC instrument the medium op- 
tical blocking filter was used, while for the MOS2 and PN de- 
tectors I ntruder etall200ll) the thin filter was chosen. All three 
X-ray cameras operated in Full Frame Mode. The data were 
processed using SAS version 5.4.1. 

A fading X-ray source was clearly detected in the IBAS 
error region. The light curve, obtained by summing the counts 
from the three EPIC cameras, is shown in Fig.|2] The afterglow 
time decay (F v (t) oc r s ) is well fit by a power-law with index 
6= 1.46 + 0.04(1 cr). An analysis of the light curves in different 
energy ranges and of the corresponding hardness ratios showed 
no evidence for spectral changes during the XMM-Newton ob- 
servation. 

Afterglow spectra from the PN and MOS cameras were ob- 
tained after standard data screening. The GRB afterglow posi- 
tion was close to a gap between two chips of the PN detec- 
tor, so the PN spectra were extracted both from a circle of 
40" radius and from a smaller circle of 20" radius, which 
did not contain the gap. Since the two extraction regions led 
to consistent results, the spectrum extracted from the larger re- 
gion, containing about 90 % of the source counts, was used for 
subsequent analysis. The background spectra were taken from 
source free regions of the same observation. The three spectra 
in the 0.4- lOkeV range, rebinned to have at least 30 counts 
per bin, are well fit by a power-law model (Fig. [3 with spec- 
tral index fix = 0.47 + 0.01 (1 cr) (where F v oc v~P"), absorption 
N# = 7.4 ± 0.9 x 1 20 cm" 2 and flux of 1 . 1 x 1 0" 1 2 erg cirT 2 s" 1 
(2- lOkeV) at 11 hours after the burst (observed flux not cor- 
rected for absorption). The column density is consistent with 
the Galactic value in this direction (-8.6 x 10 20 cm -2 ). 

From the analysis of the combined MOS and PN spectra of 
the entire observation, a 3 cr upper limit of 40 and 200 eV can 
be set on the equivalent width for narrow emission lines in the 
0.5 - 2 and 2-8 keV energy ranges, respectively. 
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Fig. 3. EPIC spectrum of the GRB 040106 afterglow and its 
best-fit absorbed power-law model. The upper data points refer 
to PN while the lower ones to MOS1 and MOS2. The residuals 
are in units of standard deviations. 



4. Discussion 

4.1 . Interpretation of the X-ray Afterglow Behaviour 

The X-ray afterglow of GRB 040106 has a remarkably hard 
fix, which is flatter than that of any other X-ray afterglow seen 
by XMM-Newton (Piro 2004). BeppoSAX measured Bx values 
in this range, but not with such high precision dFronteral20 03). 
GRB 020405 was obse rved by Chandra to have B x = 0.72 + 
0.2 1 and 6 = 1 .87 + 1 . 1 ( Mirabal e t al.l2003l) compatible, within 
the large errors, with the values measured by XMM-Newton in 
GRB 040106. 

Following the approach of IPrice et alJ d2002l) we consider 
the suitability of three afterglow models using the closure re- 
lation 6 + bBx + c — where the specific values of b and c 
required for closure depend on the model and whether the X- 
ray frequency, vx, is above or below v c , the cooling frequency 
of the electrons. The models considered are (i) a simple spher- 
ical fireball exp anding into a homogeneous medium (TSM') 
dSari et al.ll998h : (ii) a spherical fireball exp anding into a wind 
environment ('Wind') dChevalier & Lil 19991) and (iii) a spread- 
ing jet expanding into either a wind-stratified medium or the 
ISM dSari et alJl999i) . 

The combination of very hard Bx and relatively typical 
6 makes it difficult to fit the data with any of these mod- 
els (Table E). The model which comes closest to meeting the 
closure condition is the case of spherical expansion into a 
stellar-wind density profile when v c > vx- This interpreta- 
tion of the XMM-Newton d ata has already bee n considered 
in the case of GRB 040106 ddendre et alJl20()4 . The results 
of the XMM-Newton analysis presented in S ection El are con- 
sistent with those obtained by iGendre et al.l (|2Q04). Our im- 
proved constraint on 6 (+0.04 rather than +0.1) comes from 
our use of both PN and MOS camera data, improving the light 
curve statistics and resulting in a less convincing agreement 
(5 <x away from the closure condition) with the wind model. 
The wind interpretation also requires that v e is already above vx 



Table 3. Afterglow model tests using X-ray data. The con- 
dition for closure depends on whether v c > vx or v c < vx 
(Column 2). p is the power-law index of the emitting electron 
population. 



Model 




b,c 


Closure 


V 


ISM 
ISM 


> v x 
< v x 


(-3/2,0) 
(-3/2, 1/2) 


0.76 ± 0.06 
1.26 + 0.05 


1.94 + 0.02 
0.94 ± 0.02 


Wind 
Wind 


> v x 
< v x 


(-3/2,-1/2) 
(-3/2,1/2) 


0.26 ± 0.05 
1.26 + 0.05 


1.94 + 0.02 
0.94 ± 0.02 


Jet 
Jet 


> v x 
< v x 


(-2,-1) 
(-2,0) 


-0.48 ± 0.06 
0.52 ± 0.06 


1.94 + 0.02 
0.94 ± 0.02 



at the time of the X-ray measurement since no deviation from a 
power-law decay is seen in the XMM-Newton lightcurve. The 
required high cooling frequency suggests that the stellar wind 
has a low density and/or a low energy fraction for magnetic 
field s in the blast wa ve, similar to the case of GRB 020405 
dChevalier et al.l2004 . 

There are two measurements in R-band (m R _ 22.4 + 0.1 
and 2 3.7 ± 0.3) of the likely optical afterglow (M asetti et alJ 
2004), the first of which was made during the XMM-Newton 
observation, the second 24 hours later. The optical flux 
is consistent with t he extrapolation of the X-ray spectrum 
ddendre et alJl2004 . The observed temporal decay of the R- 
band light curve has a power-law index of 1.2 ± 0.4 which is 
compatible with that seen at X-ray wavelengths. Therefore, 
the cooling frequency could be below optical wavelengths or 
abov e the X-ray at the time of the XMM-Newton observa- 
tion. Gen dre et alJ (2004) argue that the cooling frequency has 
moved above the X-ray band 1 1 hours post-burst since their 
XMM-Newton analysis gives good agreement with the wind 
model prediction for this case. We have shown that there is no 
strongly compelling case in favour of the wind model from the 
X-ray data (Table E). 

In fact, the X-ray spectral index is exactly coincident with 
the 1/2 spectral index 'f ast cooling' case of the synchrotron 
shock afterglow model dSarietalJll99i IZhang & Meszaros 
2004; Piran 2000), for the limited spectral region where vx > v c 
(where v c is the synchrotron cooling frequency) and vx < v m , 
where v m represents the characteristic synchrotron frequency 
of the lower-energy limit of the emitting electron distribution. 
This limit exists so that electron energy distributions with a 
power-law index p > 2 do not diverge at low energies). The 
predicted temporal decay in this case is much flatter than that 
observed, however a more complicated geometry, such as in the 
case of an anisotropic jet, where the emission or bulk Lorentz 
factor varies with angle, leaves the time decay index largely 
unconstrained ( Meszaros et al Jl998h . 

For all cases shown in Table E] the appropriate power-law 
index for the electron distribution, p, is significantly flatter than 
the nearly universal power-law spectrum for charged particles 
accelerat ed near ultra-relativisti c shocks, which has a slope of 
2.2 - 2.3 ( Achterberg et al. 2001 ). A model of highly collimated 
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jets with flat electron spectra JDai & Chengll200ll) . although 
not meeting the closure condition imposed by the X-ray data, 
occupies the right region of the parameter space and may, with 
suitable modification, provide a good fit to the afterglow data 
in this case. 

4.2. Comparison with other INTEGRAL bursts 

Eighteen GRBs occurred within the field of view of the main 
INTEGRAL instruments up to the end of October 2004. One 
of these events, initially classified as a GRB (GRB 040903) 
is believed to be an X-ray flash (XRF) or a possible Type 
I X-r ay flare from a new transient source in the Galactic 
bulge dGotz et aljE0 04). In a second case, GRB 031203, mod- 
elling of the dust-scattered X-ray echo provided th e first evi- 
denc e for a low lum inosity, XRF source (Vau ghan et alJ2004l 
I Watson et alJl2004T) . The IBIS spectrum, however, is consis- 
tent with a single p ower-law of photon index -1.63 + 0.06 
dSazonov et alJ2004l) . typical of INTEGRAL bursts. This event 
is the only INTEGRAL burst to date for which a direct red- 
shift (z = 0.1055 + 0.0001) measurement has been made 
dProchaska et all2004 . 

A single power-law model, with photon index in the 
range -2.0 to -1.0, provides a good fit to the data for the 
vast majority of IN TEGRAL bursts in the range 20 - 200 keV 
dvon Kienlin et alJ 120031) . Only one burst detected by IBIS, 
GRB 030131, was best fit by a Band model over its whole du- 
ration, with break energy, Eo, of 70 + 20keV; photon index 
below the turnover, a, of - 1 . 4 + 0.2 and photon index above 
the turnover,/?, of -3.0 ± 1.0 dGotz et alJl2003l GRB 030131 
occurred during an INTEGRAL slew and due to telemetry re- 
strict ions only lim ited SPI data were obtained for this event 
dMoran et al. 12004. 

4.3. Constraints on redshift and energy 

A relationship between a, the y-ray spectral index below the 
spectral turnover, and redshift, derived fro m BeppoSAX y-ra y 
bursts with known redshifts, is of the form ( Amati et al. 2002): 

|or| = (2.76 ± 0.09)(1 + z )-°- 75±0 06 (1) 

This relationship reflects a dependence of a on the peak en- 
ergy, E p , due to spectral curvature, rather than a fundamental 
connection between z and a. Assuming the IBIS photon in- 
dex is a, the redshift of GRB 040106 derived from Eq. 
Z ~ 0.9+Jjj (lcr). This is a reasonable assumption since most 
B values are stee per than -2, wh ile most a values are in the 
range -2 to dPreece et aT1l2000l) . We estimate 200 keV as the 
lower limit on the break energy, since there is no curvature ev- 
ident in the INTEGRAL spectra up to this energy. In the cos- 
mological rest-frame, using the lower limit on z from above, 
Eo > 300 keV. The peak energy, E p , is given by (2 +q)x Eo and 
hence E v > 84 ke V. The relationship E p oc E m d 052 dAmati et alJ 
120021) can then be used to estimate a lower limit on the isotropic 
radiated energy, E m ^, of ~ 5 x 10 51 erg in this burst. Due to the 
weakness of this GRB, no strong constraint on the spectral lag, 



and hence on the peak luminosit y from the luminosity-lag re- 
lationship, could be determined (Norris et al. 2000J). 

The lower limit on the isotropic radiated energy is a factor 
of 10 higher than the 'mean' E y ~ 5 x 10 5l) erg determined by 
Frail et alJd2001l) when geometri c corr e ctions due to the effects 
of a jet are taken into account. iLiangl d2004l) has investigated 
the relationships between jet opening angle, 9, and the prompt 
y emission and the X-ray afterglow emission in a sample of 10 
GRBs. The 20 - 2000 keV fluence as a function of 9 is found to 
be a broken power-law, with a break at 9 = 0.1. GRB 040106 
has an estimated 20 - 2000 keV fluence of 6.8 x 10~ 8 erg cirT 2 . 
This is probably an over-estimate since it assumes the power- 
law spectrum extends up to 2 Me V, which is almost certainly 
not the case. However, it allows us to determine a lower limit 
on 9 from the y-ray data from the relation: 

S r _ 6 = 0.025 x6T 3 79 (2) 

for 9 > 0.1, where S v f, is the y-ray fluence in units of 
10 erg cm dLiangl2004l) . A lower limit of 9 > 0.22 rad is ob- 
tained from this relation. Hence, a lower limit of 1.2 x 10 50 erg 
on the energy radiated in y-rays is derived. Uncertainties in 
the cosmological correction which should be applied, and in 
the true spectral shape, probably give rise to an uncertainty of 
order 2 in the y-ray fluence, which is ~ 50% on 9. 

The 9 value pr edicted from the X-ray afterglow decay slope 
ofa dLiangl2004l) according to the relation: 

6 = (1.6 ± 0.07) -(2.27 ±0.54) X 9 (3) 

is ~ 5 x 10~ 2 , considerably smaller than the lower limit derived 
from the y-ray data. If this correction is applied to the isotropic 
radiated energy estimate, then a lower limit of 6 x 10 48 erg on 
the energy radiated in y-rays is obtained. 

A third relationship, between the rest-frame X-ray after- 
glow flux 10 hours after the burst and 9, is given by: 

F Xt _ 13 = 0.47 x r L8±a4 (4) 

where F x -i3 is the rest-frame 2 - lOkeV flux derived from the 
observed flux in units of 10~ 13 ergcrrT 2 s _1 according to: 

F x = ^(l+zf^ (5) 
1 +z 

Using z ~ 0.9 derived above, a 9 value of 0.24 rad is ob- 
tained, consistent with the estimate from the y-ray data. 

From the burst morphology view-point, GRB 040106 has 
the characteristics of a low-luminos ity, long lag burst (few, 
well-separated peaks) (Norris 2002). Its location in super- 
galactic coordinates (161°, - 22°) is consis tent with that found 
in the subset of such bursts (Norris l2002l) . but there is no in- 
dependent evidence to suggest that this event is anything other 
than a 'standard' cosmological burst. The pulse properties and 
time i ntervals between pulses are related to T90 dMcBreen et alJ 
2002). The time interval, AT, and properties of the two pulses 
of GRB 040106, including rise time, fall time and FWHM fit 
well w ith the expected trends from previous analyses of BATSE 
bursts (McBree n et al.l2002l) . 
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Fig. 4. Extrapolated 2-10keV fluxes from SPI (circles) and 
IBIS (squares) for each of the two pulses. The solid line is the 
extrapolation of the X-ray data while the dashed lines repre- 
sent the range of decay slopes based on the power-law decay 
of the 2 nd pulse. 



4.4. GRB-Afterglow Connection 

We examined the possibility that the second pulse of the GRB 
is the onset of the afterglow. The temporal decay of the 2 nd 
pulse is consistent with a power-law (with large errors due to 
the small number of bins and their low counts) which may in- 
dicate the presence of a high-energy afterglow, due to external 
shocks, during the burst itself. 

The IBIS and SPI derived fluxes in the two pulses were 
extrapolated to X-ray energies, assuming a single power-law 
spectrum, to estimate the 2-10keV flux during the GRB. 
Clearly this is valid only if the GRB has neither an X-ray ex- 
cess nor is X-ray rich. It is unlikely that this is an X-ray rich 
GRB with a peak energy at or below the low end of the SPI and 
IBIS detector sensitivity (i.e. ~ 20keV) since the photon index 
would then correspond to an unusually har d value for the high - 
energy index above the spectral turnover (Pree ce et alJ l2000). 
It is more likely that the weakness of this GRB washes out ev- 
idence for a spectral break at more typical energies of a few 
hundred keV. 

The extrapolated fluxes from SPI and IBIS are shown as 
symbols in Fig. |4] Folding in the range of possible values of 
the SPI and IBIS extrapolated fluxes, assuming they follow the 
power-law temporal decay of the 2 nd GRB pulse, these extrap- 
olated values fall marginally below the backward extrapola- 
tion of the late-time X-ray afterglow data obtained by XMM- 
Newton (Fig. [4}. Assuming the decay slope of the 2 nd pulse is 
indicative of the onset of a high-energy afterglow, the change 
in temporal decay occurs roughly between 50 and 5 000 s after 
the onset of the GRB and may be evidence of the passage of the 
cooling break through the X-ray band in that time window. The 
fact that the y-ray photon index and X-ray spectral index are 
within 2 <x of each other suggests that v c < vx during the X-ray 
observations. This interpretation hinges on the assumption that 
the decay of the 2 nd GRB pulse is the onset of the afterglow. 



There is no evidence in the IBIS/ISGRI light curve for 
soft extended or delayed emission such as that observed by, 
for ex ample, SIGMA/GRANAT in GR B 920723 feurenin et all 
1999) or by HETE-II in GRB 02121 1 dCrew et all2003l) . 

5. Conclusion 

INTEGRAL'S capabilities for GRB studies, from the rapid and 
accurate localisation provided by IBAS, to the good timing and 
spectral resolution provided by SPI and IBIS, are well illus- 
trated in the case of GRB 040106. The y-ray measurements 
have provided a constraint on the redshift (z ~ 0.9^ ^) and en- 
ergy (lower limit on the isotropic radiated energy of ~ 5 x 10 51 
erg) for this burst, in the absence of optical spectroscopy. 
Furthermore, previously identified correlations between y-ray 
properties and X-ray afterglows have been used to estimate 
the jet opening angle (~ 0.22 rad) and hence the geometry- 
corrected radiated energy (1.2 x 10 50 erg) in this burst. 

The combination of very hard energy index and relatively 
typical temporal decay makes it difficult to fit the X-ray data 
with afterglow models of a simple spherical fireball or a jet, ex- 
panding into either the ISM or a wind environment. In particu- 
lar, we can rule out the case of isotrop ic expansion into a wind 
environment (v c > vx), suggested by iGendre et alJ d2004l) . at 
5 <x confidence. In all scenarios, the power-law index of the 
emitting electrons is flatter than 2 and models involving flat 
electron spectra may be more successful fitting the X-ray af- 
terglow data of this burst. 

The additional science provided by the combination of y— 
ray and X-ray early afterglow measurements highlight the im- 
portance of broadband high-energy coverage from the GRB 
through to the first days post-burst. We anticipate the signifi- 
cant progress in this area that will be made with the launch of 
NASA's Swift mission. 

Acknowledgements. We thank the anonymous referee for his/her use- 
ful comments which helped to improve this paper. 
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